x_max = 2*pi;
N = 16;
mu = -0.5;
max_iter = -N/mu;
numfine = 3;
FTCS_errmax = zeros(1,numfine+1);
for i = 1:numfine+1
    u = initial_data(N,x_max);
    u = FTCS(u,N,mu,max_iter);
    exact_u = exact_solution(N,mu,x_max,max_iter);
    FTCS_errmax(i) = err(u,exact_u,-1,x_max);
    N = N*2;
    max_iter = max_iter*2;
end
printOrder(FTCS_errmax);